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We study the phase space dynamics of multi-dimensional symplectic maps, 
using the method of the Generalized Alignment Index (GALI). In particular, 
we investigate the behavior of the GALI for a system of N = 3 coupled standard 
maps and show that it provides an efficient criterion for rapidly distinguishing 
between regular and chaotic motion. 
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1. Introduction 

The distinction between regular and chaotic motion in conservative dy- 
namical systems is fundamental in many areas of applied sciences. This 
distinction is particularly difficult in systems with many degrees of freedom, 
basically because it is not feasible to visualize their phase space. Thus, we 
need fast and accurate tools to obtain information about the chaotic vs. 
regular nature of the orbits of such systems and characterize efficiently 
large domains in their phase space as ordered, chaotic, or "sticky" (which 
lie between order and chaos). 

In this paper we focus our attention on the method of the General- 
ized ALignment Index (GALI), which was recently introduced and ap- 
plied successfully for the distinction between regular and chaotic motion in 
Hamiltonian systems. 1 The GALI method is a generalization of the Smaller 
Alignment Index (SALI) technique of chaos detection. 2-4 We present some 
preliminary results of the application of GALIs on the dynamical study 
of symplectic maps, considering in particular the case of a 6-dimcnsional 
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(6D) system of three coupled standard maps. 5 It is important to note that 
maps of this type have been extensively studied in connection with the 
problem of the stability of hadron beams in high energy accelerators, see 6 
and references therein. 

Our numerical results reported here verify the theoretically predicted 
behavior of GALIs obtained in 1 for Hamiltonian systems. In addition we 
study in more detail the behavior of chaotic orbits which visit different 
regions of chaoticity in the phase space of our system. 

2. Definition and behavior of GALI 

Let us first briefly recall the definition of GALI and its behavior for regular 
and chaotic motion, adjusting the results obtained in 1 to symplectic maps. 
Considering a 2A-dimcnsional map, we follow the evolution of an orbit 
(using the equations of the map) together with k initially linearly indepen- 
dent deviation vectors of this orbit T?i, 1?2, with 2 < k < 2N (using 
the equations of the tangent map). The Generalized ALignmcnt Index of 
order k is defined as the norm of the wedge or exterior product of the k 
unit deviation vectors: 

GALI k (n) HI i>i(n) A v 2 {n) A ... A v k (n) || (1) 

and corresponds to the volume of the generalized parallelepiped, whose 
edges are these k vectors. We note that the hat (') over a vector denotes 
that it is of unit magnitude and that n is the discrete time. 

In the case of a chaotic orbit all deviation vectors tend to become linearly 
dependent, aligning in the direction of the eigenvector which corresponds 
to the maximal Lyapunov exponent and GALI^ tends to zero exponentially 
following the law: 1 

GALI k {n) oc e -[(i-^)+(i-<*)+-»+(*i-«'*)]" ) (2) 

where <7i , . . . , a k are approximations of the first k largest Lyapunov expo- 
nents. In the case of regular motion on the other hand, all deviation vectors 
tend to fall on the A-dimensional tangent space of the torus on which the 
motion lies. Thus, if we start with k < N general deviation vectors they 
will remain linearly independent on the A-dimensional tangent space of 
the torus, since there is no particular reason for them to become aligned. 
As a consequence GALL remains practically constant for k < N. On the 
other hand, GALL; tends to zero for k > N, since some deviation vectors 
will eventually become linearly dependent, following a power law which de- 
pends on the dimensionality of the torus on which the motion lies and on 
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the number m (m < N and m < k) of deviation vectors initially tangent 
to this torus. So, the behavior of GALI^ for regular orbits is given by 1,7 

{constant if 2 < k < N 
n . i{k } N) _ m if N < k < 2N and < m < k - N . (3) 
if N < k < 2N and m > k - N 

3. Dynamical study of a 6D standard map 

As a model for our study we consider the 6D map: 

X ^ — X ]_ ~\~ x^ 

x' 2 = x 2 + Jjp sin(27rxi) - ^r{sin[27r(x 5 — .ti)] + sin[27r(a; 3 — a;i)]} 
x' 3 = x 3 + x'^ 

x' A = x± + -J! sin(27ra;3) - J;:{sin[27r(xi - x 3 )\ + sin[27r(x 5 - x 3 )]} 
x r-> = x 5 + x e 

x' 6 = x e + £ sin(27rx 5 ) - ^ F {sin[27r(xi - x 5 )] + sin[27r(x 3 - x 5 )}} 

which consists of three coupled standard maps 5 and is a typical nonlinear 
system, in which regions of chaotic and quasi-periodic dynamics are found 
to coexist. Note that each coordinate is given modulo 1 and that in our 
study we fix the parameters of the map (4) to K = 3 and /5 = 0.1. 

In order to verify numerically the validity of equations (2) and (3), 
we shall consider two typical orbits of map (4), a chaotic one with initial 
condition x\ = x$ = X5 = 0.8, x 2 = 0.05, X4 = 0.21, x$ = 0.01 (orbit CI) 
and a regular one with initial condition x\ = x% = X5 = 0.55, x 2 = 0.05, 
x 4 = 0.01, xq = (orbit Rl). In figure 1 we see the evolution of GALIfc, 
k = 2, . . . , 6, for these two orbits. 

It is well-known that in the case of symplcctic maps the Lyapunov 
exponents are ordered in pairs of opposite signs. 8 Thus, for a chaotic orbit 
of the 6D map (4) we have o~\ = —ctq, a 2 = —0-5, as = — 04 with a\ > a 2 > 
03 > 0. So for the evolution of GALIfc, equation (2) gives 

GALI 2 (n) oc e-^ 1 "" 72 )™, GALI 3 (n) oc e -(^i-^-<r 3 )n ^ 
GALI 4 (n) oc e-^ 1 -^)™, GALI 5 (n) oc e - 4CTl ™, GALI 6 (n) oc e - 6<Tl ™. 

(5) 

The positive Lyapunov exponents of the chaotic orbit CI were found to 
be <Ti w 0.70, (72 w 0.57, (73 w 0.32. From the results of figure la) we see 
that the functions of equation (5) for a\ = 0.70, a 2 = 0.57, cr 3 = 0.32 
approximate quite accurately the computed values of GALIs. 

For the regular orbit Rl we first considered the general case where no 
initial deviation vector is tangent to the torus where the orbit lies. Thus, 
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Fig. 1. The evolution of GALI^, k = 2, . . . ,6, with respect to the number of iteration n 
for a) the chaotic orbit CI and b) the regular orbit Rl. The plotted lines correspond to 
functions proportional to e -( CT i- CT 2)«, e -(2<T 1 - CT2 -<T 3 )n ) e -(3 CT1 - CT2 )n _ e -4o-in ; e -6a in 

for <t i = 0.70, fJ2 = 0.57, 0-3 = 0.32 in a) and proportional to n~ 2 , n -4 , ra~ 6 in b). 

for the behavior of GALI^, k = 2, . . . , 6, equation (3) yields for m = 

GALI 2 (n) oc constant, GALI 3 (n) oc constant, GALL^n) oc 4y, . . 
GALI 5 (n)oc^, GALI 6 (n)cx^. 

From the results of figure lb) we see that the approximations appearing in 
(6) describe very well the evolution of GALIs. 

In order to verify the validity of equation (3) for 1 < m < 3, in the case 
of regular motion, we evolve orbit Rl and three random initial deviation 
vectors for a large number of iterations (in our case for 5 x 10 7 iterations), 
in order for the three deviation vectors to fall on the tangent space of the 
torus. Considering the current coordinates of the orbit as initial conditions 
and using m = 1 or m = 2 or m = 3 of these vectors (that lie on the tangent 
space of the torus) as initial deviation vectors we start the computation of 
GALIs' evolution. We note that the rest 6 — m initial deviation vectors 
needed for our computation are randomly generated so that they do not 
lie on the tangent space of the torus. The results of these calculations are 
presented in figure 2, where the evolution of GALI^, k = 2, ... ,6, for dif- 
ferent values of m is plotted. Figure 2 clearly illustrates that equation (3) 
describes accurately the behavior of GALIs for regular motion also in the 
case where some of the initial deviation vectors are chosen in the tangent 
space of the torus. 
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Let us now consider the case of a chaotic orbit which visits different 
regions of chaoticity in the phase space of the map. The orbit with initial 
conditions x\ = x 3 = x 5 = 0.55, x 2 = 0.05, X4 = 0.21, x e = 0.0 (orbit C2) 
exhibits this behavior as can be seen from the projections of its first 1000 
successive consequents on different 2-dimensional planes plotted in figure 
3. 




Fig. 3. Projections on the planes a) (xi, 2:2), b) (2:3, X4) and c) (x5,xe) of the first 1000 
successive points of the chaotic orbit C2. 



The projections look erratic, indicating that the orbit is chaotic. How- 
ever, we also observe in all three projections that the C2 stays 'trapped' 
for many iterations in two oval-shaped regions and eventually escapes en- 
tering the big chaotic sea around these regions. This behavior is also de- 
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picted in the evolution of the Lyapunov exponents of the orbit (see figure 
4a). The three positive Lyapunov exponents are seen to fluctuate around 
(Ji « 0.033, cr 2 ~ 0.02, cr 3 w 0.005 for about 1000 iterations exhibiting 
'jumps' to higher values when the orbit enters the big chaotic sea, stabiliz- 
ing around a x « 0.793, cr 2 ~ 0.624, cr 3 « 0.365. 




Log(n) n 



Fig. 4. a) The evolution of the three positive Lyapunov exponents of the chaotic orbit 
C2. The evolution of GALI fc , k = 2, . . . , 6, with respect to the number of iteration 
n for the same orbit when we use as initial condition of the orbit its coordinates at 

b) n = and at c) n = 10 iterations. The plotted lines in b) and c) correspond to 
functions proportional to e -( CT i- CT 2)". j e -(2a 1 -a 2 -a 3 )n ^ e -(3c 1 -c 2 )n ^ e -4<Tin 5 e -6<rm 

for (Ti = 0.033, (72 = 0.02, ct 3 = 0.005 in b) and a 1 = 0.793, a 2 = 0.624, cr 3 = 0.365 in 

c) . 

Let us now study how the GALIs are influenced by the fact that orbit 
C2 visits two different regions of chaoticity characterized by different values 
of Lyapunov exponents. Since C2 is a chaotic orbit, its GALIs should tend 
exponentially to zero following the laws of equation (5). Thus, starting 
the computation of the GALI&, k = 2, . . . , 6, from an initial point of C2 
located in the first chaotic sea, we see that the slopes of the exponential 
decay of GALIs are well described by equation (5) using for <ti, ct 2 , the 
approximate values of the Lyapunov exponents of the small chaotic region 
(figure 4b). On the other hand, using as initial condition for this chaotic 
orbit its coordinates after 10 6 iterations, when the orbit has escaped in the 
second chaotic region, the evolution of GALIs is again well approximated 
by equation (5) but this time for u\ = 0.793, CT2 = 0.624, a 3 = 0.365, which 
are the approximations of the Lyapunov exponents of the big chaotic sea 
(figure 4c). Thus, we see that in this case also the Oi, i = 1, . . . , k which 
appear in equation (5) are good approximations of the first k Lyapunov 
exponents of the large chaotic region in which the orbit eventually wanders 
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after about 10 3 iterations. 

4. Conclusions 

In this paper we verified the theoretically predicted behavior of the Gen- 
eralized Alignment Index (GALI) by considering some particular regular 
and chaotic orbits of a 27V-dimensional symplectic map with N — 3. In 
particular, we showed numerically that all GALIfe, 2 < k < 2N, tend to 
zero exponentially for chaotic orbits, while for regular orbits they remain 
different from zero for 2 < k < N and tend to zero, following particular 
power laws, for N < k < 2N. Thus, by using GALI& with sufficiently large 
k, one can infer quickly the nature of the dynamics much faster than it is 
possible by using other methods. 

Also, the study of chaotic orbits which visit different regions of chaotic- 
ity in the phase space of the system, provides further evidence that the 
exponents of the exponential decay of GALIs are related to the local values 
of Lyapunov exponents. Thus, we have shown that the different behaviors 
of the GALIs for regular and chaotic orbits can be used for the fast and 
accurate identification of regions of chaoticity and regularity in the phase 
space of symplectic maps. We plan to investigate this further in a future 
publication concerning maps with N >> 3, which describe arrays of con- 
servative nonlinear oscillators. 
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